Multi-electron SEFs for nuclear reactions involved in advanced 

stages of stellar evolution 

Theodore E. Liolios ^'^'^ * 
^European Center for Theoretical Studies in Nuclear Physics and Related Areas 
Villa Tamhosi, 1-38050 Villazzano (Trento), Italy 
University of Thessaloniki, Department of Theoretical Physics 
Thessaloniki 54006, Greece 
3 Hellenic War College, EST 903, Creece 



Abstract 



Multi-electron screening effects encountered in laboratory astrophysical re- 
actions are investigated by considering the reactants Thomas-Fermi atoms. 
By means of that model , previous studies are extended to derive the corre- 
sponding screening enhancement factor (SEF), so that it takes into account 
ionization, thermal, exchange and relativistic effects. The present study, by 
imposing a very satisfactory constraint on the possible values of the screening 
energies and the respective SEFs, corrects the current (and the future) exper- 
imental values of the astrophysical factors associated with nuclear reactions 
involved in advanced stages of stellar evolution. 
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I. INTRODUCTION 



The astrophysical factors for the nuclear reactions of the CNO bi-cycle (solar center, 
upper main-sequence stars etc.) and of other advanced nuclear burning stages, such as 
those ignited in classical novae, have been obtained by performing measurements well above 
the Gamow peak, ( fl]] and references therein) and extrapolating to lower energies thus 
committing a significant error in certain cases. After the accomplishments of the LUNA 
collaboration in probing deep into the Gamow-peak region of the most important break- 
up reaction He^ {He^, 2p) He^ the nuclear astrophysics community hopes that similar low- 
energy experiments will follow for astrophysical nuclear reactions involving heavier nuclei. 
However, the screening effect, which has been studied extensively ( 0,^] and references 
therein) regarding nuclear reaction of the proton-proton chain, is particularly accentuated 
when heavier ions are considered due to their multi-electron nature. It is the purpose of 
this paper to investigate the enhancement of the astrophysical factor that appears (or will 
appear) when astrophysical energies are attained in nuclear reaction experiments. 

In a recent paper [Q] the Thomas-Fermi (TF) model was used in order to derive a lower 
(sudden) limit (SL) and an upper (adiabatic) limit (AL), thus imposing a constraint on the 
possible values of the screening energies. In the present work we will improve that study by 
tightening that constraint while studying, at the same time, ionization, relativistic, thermal 
and exchange effects. 

The TF model is one of the most thorougly studied models of atomic physics (for 
a detailed study see Ref. [T^). In brief, if n (r) is the electron number density around 
the screened nucleus then the self- consistent potential $ (r) , which results from the com- 
bined field (nucleus plus electrons), is obtained by solving the Poisson equation (r) = 
—Anp (r) where p (r) = —en (r) is the charge density of the electron cloud. Using quantum 
statistical arguments |jTO| we can show that for the TF atom at each point r the fermi 
energy pj must be 

^ = e [$ (r) - $ (ro)] (1) 
where ro is the ionic radius which gives the boundary conditions 

nr.)-^,(f) (2) 

ro \dr J ro 

and the degree of ionization is a function of the number of electrons (A^) and protons (Z) , 
respectively, so that q = 1 — N/Z. 

Following some well cited transformations |jlO| the self-consistent potential can be writ- 
ten: 

$ (r) = $ (ro) + —0 (3) 
r 

Setting r = ax, the universal function (x) is given by the dimensionless TF equation: 
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(4) 

to be solved with the boundary conditions (f) (0) = 1, (j)' (xq) = —q/xo- Note that the scahng 
parameter a has its usual value a = 0.88534:Z~^^'^aH, with representing the Bohr radius. 



II. IONIZATION, EXCHANGE AND THERMAL EFFECTS 

We will exploit the advantages of the TF theory in order to derive the effects of ion- 
ization, electron-electron (exchange) interactions and temperature on multi-electron SEFs. 
The adiabatic limit (AL) SEF of Ref. |0 for neutral atoms seems to be well justified and 
it is the largest available. However, the sudden limit (SL) SEF of the same work needs 
improvement as it has been derived using Tietz's approximation [jll| which is actually a 
rough model that can only approximate the self-consistent potential of neutral atoms at 
average distances from the nucleus, disregarding all relativistic, thermal and ionization 
effects. As far as astrophysical reaction are concerned, screening becomes important very 
close to the nucleus and therefore a more appropriate approximation is needed, which could 
also account for ionization effects. To solve that problem we resort to Baker's small 
X expansion where the expansion coefficients are themselves functions of the initial slope 
S {q) = (0) • The first few terms are given by: 

<Pb{x) = 1 + S {q) X + ^x3/2 + Is (q) x'/' + ^x' + ^S' (q) x' + ^5 (g) x^ + ... (5) 

3 5 3 7U 15 

There exists a unique value 5* (0) = —1.588 which gives xq = cxd and (xq) = correspond- 
ing to the neutral TF atom (g = 0). Every other value S {q) < S (0) gives finite values 
for Xq corresponding to positive TF ions (0 < g < 1) . By selecting various initial slopes we 
obtain relations of the form: S = S (q) and Xq = Xq (g) (see Figures 1 and 2) 

Hence, the screened Coulomb potential around the target nucleus will now be given by: 

'^W = ^A + — - (6) 



axo (g) r \a 

where from now on Zi and Z2 represent the atomic numbers of the target nucleus and the 
projectile, respectively 



An alternative approach would be to adopt Moliere's approximation [jT3[ for the solution 
ofEq. d): 

3 

{x) = ^Bi exp (-pix) (7) 
1=1 

with the fitting parameters 

Bi = 0.1, B2 = 0.55, 53 = 0.35, pi = 6.0, P2 = 1-2, Ps = 0.3 (8) 



3 



However, Moliere's function is less accurate than Baker's one since the former un- 
derestimates the initial slope. For example, for neutral atoms (pj^^ (0) = —1.365, while 
(0) = -1.588 

When positive ions are considered the total binding energy |Tg Et^^p = -20.93 Z'^/^ eV 
used in Ref. [01 for neutral TF atoms must be modified to account for ionization. In such 



a case 



where 



F{q) 
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S{q) + 



Xq (q) 



(9) 



(10) 



We can now derive the screening enhancement factor for astrophysical nuclear reactions. 
The energy dependent penetration factor P (E) multiplied by the astrophysical factor S (E) 
in the s-wave cross section formula 



is given by the WKB method: 

(E) = exp 



E 



2^ r^E) 



h 



(r) - Edr 



(12) 



where (SC) stands for screening and the classical turning point Vc is given by 

Vsc (rj = E (13) 

while the nuclear radius r„ is a function of the mass numbers ^41^2 of the reactants, so that 
r„ = 1.4(<VA^/=^)/m. 

At astrophysical energies in the laboratory the potential energy is found to be shifted 
by a constant screening energy Ug which is added to the relative energy of the collision. For 
the no-screening (NOS) nucleus-nucleus channel the calculation is easy leading to [0 



pNOS _ 



-27m 1 



4 /r„\i/2 



+ 



3tt \r, 



3/2 



+ 



where n is the Sommerfeld parameter and Ug = 0. The general SEF / (E) is defined as 



(14) 



fiE) 



pNOS ^E) 



(15) 



which for weak screening (WES) shifts Ue « E can be written ||T7 

U, 



fWES ^ , 



E 



(16) 



4 



In fact all measurements in the laboratory should be corrected by means of the above SEF 
so that the bare-nucleus astrophysical factor is not overestimated. Following the treatment 
of Ref. which was used for neutral TP atoms, we can derive an approximate analytic 
formula for the TF screening shifts U^^J"'^^^ and the corresponding SEFs f^F'^^^ i-^) 
the SL and the AL respectively, where the degree of ionization is also taken into account. 
Combining Eqs. (16) , (12) ,and (6) we obtain the following formulas: 
Sudden Limit 



Siq) + 



e 

an 



(17) 



/#(E)^exp 



0.482Z[/'Z|Ai/2 



E 



3/2 
(keV) 



Xo (q) 



(18) 



Adiabatic Limit 



AL 
TF 



F (qu) {Z, + Z,y/^ - Fi (gi) Zi^' - F, (q^) Z^ 



7/3 



.7/3" 



(19) 



/#(E)^exp 



15.64^1^2 F(gi2)(Zi + Z2 



,7/3 



Fliqi 



zT - 



F2 {qi 



zr 



^1/2 ■ 



^3/2 



(20) 



where the subscripts (1,2,12) of q in the adiabatic limit formulas take into account the 
different degrees of ionization of the reactants and the combined nuclear molecule. 

The present sudden limit formulas are more accurate than the ones obtained by means of 
Tietz's approximation, therefore improving the constraint on the possible range of screening 
energies. For example the lower limit for neutral atoms given by Eq. (^) is 40% larger 
than the corresponding one of Ref. , thus greatly narrowing the gap between the sudden 
and the adiabatic limit SEFs. 

Further simplifications are possible by noting that, for small ionization numbers g, the 
second additive term q^ /xq (g ) in Eq. (|10D is much smaller than the initial slope 5* (g) and 
can therefore be neglected. If we further assume proton induced nuclear fusion reactions 
where the ionization state before and after the collision is the same then the AL formulas 
can be written as 



TjAL _ — 



12 

7 



1/3 



+ 1)'/' - zi 



S{q)^ 



q^ 

a;o (g) 



e 

an 



(21) 



/#(E)^exp 



\S (g) + gVo^o (g)] 

5^3/2 



Z1Z2 (Zl + 1)"/=^ - 



.7/3 



^1/2 



:22) 



Note that if we set g = in all the above SL and AL formulas we recover the respective 
ones for neutral targets . However, the jj^p (E) SEF given by Eq. (|1^) must be used 
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with caution. It has been derived by truncating Baker's series but at very low energies and 
moderately heavy nuclei that approximation begins to break down. Its actual validity will 
be investigated along with the relativistic effects in Sec. III. 

We can assess the quality of the constraint imposed by our model on the possible values 
of the screening energies and the associated SEFs by studying Fig. 3 and Fig. 4. In the 
former, the SL and AL screening energy U^p'^^'' for the most important astrophysical 
nuclear reactions of the CNO bi-cycle are plotted with respect to the degree of ionization 
of the TF positive ionic target . Along the whole profile of the q values the upper {AL) 
and lower {SL) screening energy limits differ by roughly 100 eV (or less) which is a very 
satisfactory constraint for most practical purposes. 

On the other hand in Fig. 4 we can observe the SL/AL SEFs for the most important 
astrophysical nuclear reaction of the CNO bi-cycle, namely A^^^ (p, 7) O^^, with respect to 
the center-of-mass energy E in the region of the solar Gamow peaks {Egp = 27.2 keV) for 
various degrees of ionization. Naturally, the more ionized the atomic target the smaller the 
screening energy (see also Fig. 3) and, consequently, the smaller the respective SEF at a 
particular energy. The quality of the SL/AL constraint on the SEF is very good as has 
already been indicated by the corresponding screening energy limits (see Fig. 3.) 

Before we proceed to study relativistic effects, we should consider the influence of ther- 
mal and exchange effects. It should be emphasized that we have defined the sudden limit 
for a cold electron gas (T = 0) disregarding all exchange effects which, on the contrary, 
have been taken into account in the derivation of the total TF binding energy Ei^p {q) . Of 
course, the latter also includes the total potential energy Ej!'p^ between the nucleus and the 
electrons and the total kinetic energy Kpp of the bound electrons. Actually, the exchange 
energy Elfphs expected to lower the relative collision energy as recently shown in a relevant 
study . For a neutral atom the above components of the total binding energy Eifp {q) 



can be written W8 



= +20.99Z^/^ E^TF^ = +7.00Z^/^ E^f = -48.99^^/^ (23) 

where we observe that exchange effects play a less important screening role than the kinetic 
energy ones. 

It is now obvious that the adiabatic limit SEF can be written as 

/# {E) = f^L {E) X f^t., {E) X f^i.^, {E) (24) 

where the first two factors of the right-hand side decrease rather than increase the cross 
section of the astrophysical fusion reaction. 

The above factorization of multi-electron SEFs allows some qualitative investigation 
of the phenomenon. Obviously, according to the third factor, by increasing the binding 
energy of the electrons we can enhance the cross section of the reaction . That is for the 
same number of electrons, the larger the atomic number the larger the cross section at a 
particular collision energy. However, by increasing the atomic number the kinetic energy 
of electrons is also increased and according to the first factor that will cause a decrease in 
the cross section. On the other hand exchange effects also increase with the atomic number 
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reducing further the cross section. Note that the dependence of the SEF on the atomic 
number is more complex than the simple monotonous increase indicated by Eq. (^) and 
Eq. (^3|) . The TP model, being a statistical one, gives a net dependence of the total energy 
on the atomic number proportional to Z'^l'^ . Unfortunately this conceals all shell effects, 
which are very important, as each shell carries a different screening capacity. For example, 
the K-shell electrons are expected to play a much more important role than those lying in 
the L-shell and so on. In fact the sudden limit SEFs which were derived Q for the first and 
second excited states in hydrogen-like atoms, showed that a single electron of the K-shell 
increases the relative collision energy four times more than a single electron of the L-shell. 
As the screening energy enters the SEF exponentially (Eq. ([16|) ) that translates into a 
K-shell SEF being equal to the L-shell SEF raised to the fourth power. Higher shells can be 
shown to play an even minor role. Notably, the present TF approach for the sudden limit 
disregards the fast response of the electron cloud due to the combined nuclear molecule p , 
thus obtaining an even more conservative value for the SL screening energy. 

In order to study the relative importance of shell effects we can consider the screening 
energy generated by the collision of a bare Carbon nucleus (Z = 6) with an atom of the 
same nuclear charge with only two electrons in its closed K-shell. Disregarding the fast 
response of the two-electron cloud, so that we are consistent with the present TF model, 
the Hartree-Fock (HF) screening energy is P| 



TjSL 



5 _ . e2 



-2Z {Z - 5/16) + -{Z - 5/16) 



an 



(25) 



with a numerical value of = 1760 eV. 

On the other hand the Thomas- Fermi SL model (Eq.(|l^) ) for g = gives Uj'p = 
3191 If we also take into account the exchange effects disregarded by Eq. (plTj) the TF 
screening energy should be smaller. In any case, according to our models, more than 50% 
of the six-electron TF screening energy originates from the two K-shell electrons, which 
emphasizes the importance of shell effects. 

The TF formulas derived in this paper show a direct dependence of multi-electron SEFs 
on the ionization state (g) of the reacting atoms. The ionization number q is of course 
a function of temperature (see for example the Saha equation []I6[) so that at very high 
temperatures the atom is completely ionized (e.g. stellar plasma). Thus, the cold electron 
gas assumption is no longer valid when the temperature rises. For thermal energies of a 
few eV , corresponding to a few thousand degrees Kelvin, Eq. (H) should be replaced by 
the Marshak-Bethe equation ||19||: 



'MB 



(x,T) 



X 



1/2 



1+C(T) 



X 



(26) 



where the dimensionless temperature-dependent parameter is ( (T) = lO'^^Z^ , with 

T measured in degrees Kelvin. The solution of that equation can be found by a perturbation 
method |0 and reads: 



(f)MB (X, T) = 00 (x) + C (T) 01 (x) 



(27) 
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where the unperturbed term 0o i^) is the solution of Eq. and the trial function (pi (x) 
is approximated by the truncated series 0i (x) ~ (4/35) x''/^ — (4/63) S (0) x^/^ + ... Conse- 
quently, in our study there is now a temperature dependent screening correction Ue (T) to 
be taken into account. After some algebra that term for small distances can be written 

UJT) = ±((T)^(Lf (28) 



35 a \aj 



Actually, in obtaining Eq. (|T^ we disregarded radially dependent terms arising from 
Baker's expansion since even the largest of those neglected terms, given by 

^SL, ^ 4ZiZ2e2 /r\V2 



U^'^{r) = --^^{-\ (29) 
3 a \aj 

is much smaller than the constant screening shift throughout the tunnelling region. By 
comparing Eq. (^) and Eq. (p9|) it is obvious that for temperatures of several thousand 
degrees Kelvin (so that kT « e$ (ro) ) the multi-electron screening enhancement of 
astrophysical reactions is the same as in a cold electron gas. 



III. RELATIVISTIC EFFECTS 



We can now proceed to include relativistic corrections to screening energies and the 
respective SEFs, which are expected to become important when studying the screening 
effects of lower shell electrons on a heavy atomic nucleus. For example a K-shell electron 
of an Oxygen atom has a Bohr velocity of Ue = 0.058c, which indicates that relativity 
may further modify the screening effect and therefore warrants further investigation. For 
simplicity, in the study that follows we will focus on neutral targets, as any degree of 
ionization will simply diminish the associated relativistic corrections. 

The sudden limit can be obtained by using the Dirac-Hartree-Fock-Slater (DHFS) cal- 
culations of Salvat et al 
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which provide us with screening functions accounting for 
relativistic effects. In that work Moliere's approximation was used and the corresponding 
Bi (Z) , f3i (Z) parameters were tabulated for Z = 1 — 92 . Therefore, when relativistic 
effects are considered, the SL screening shift can be written as 



TjSL 
^ DHFS 



-4 



/ 2 



V97r2 



1/3 



Y,B,{Zi)l3:,{Zi 



.1=1 



e 

an 



(30) 



Unfortunately, Eq. (|30D is an very rough approximation. Moliere's approach was previously 
shown to underestimate the non-relativistic (neutral atom) SL screening energy by 14%. 
In our study we have used the relativistic and non-relativistic screened Moliere potentials 
given in Ref. |]l3l and Ref. [^], respectively. Although they both satisfy the condition that 
relativistic effects result in more bound systems none of them is reliable close to the nucleus. 
For example Eq. (|30|) for the reaction C^^ {p, 7) N^"^ yields U^^pg = 406e\^ which is even 
lower than the non-relativistic SL screening energy given by Eq. (0). Therefore, the DHFS 
screening energy turns out to be inappropriate as a lower bound of the relativistic screening 
energy. 
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The most accurate way to calculate the screened Coulomb potential would be to use 
the Vallarta-Rosen relativistic TF equation p2[ 



X 



1 + A 



4>rel {x) 



X 



3/2 



(31) 



with A = 3 X 10 ^Zi^^ . The V-R potential has been approximated with respect to the 
distance Xc = A + S* (0) A^, by a two-branch function . Namely, for small x and x > Xc 



€e! i^) = E (^)] - 1-206A3/2 (1 + 0.553 In Z) 



(32) 



n=0 



where 0o (x) is the non-relativistic Baker's expansion given by Eq. (|^) and 

01 (x) = -6x^/2 + 11.733a; + 5S (0) x^/^ + ... (33) 



^2 {x) = ^x'/'-^S{0)x'/' + ... 



(34) 



The only significant terms added to Baker's expansion through the sum of Eq. (32) are 
those raised to powers n = 0, 1. The remaining terms (n > 1) are negligible throughout 
the tunneling region as they are multiplied by the very small number A" and therefore 
vanish at distances smaller than the classical turning point (see also previous section and 
the Appendix in Ref. 0. 

For distances shorter than Xc (but larger than the minimum internuclear distance given 
by Xn = rn/a) the approximation is 



3 

(t)lli (x) = l + sx + -A^/2 [{xn -x)+x\n (x/xn)] + A^/2 



In 



X 



X'n 



where s is given by 

s = S{0)- X^'^ hiZ- 0.1388A^/2 _^ 11.733^ + 5.1575 (0) X^l^ + .. 



+ ... (35) 



(36) 



Instead of seeking an approximate analytic formula for the relativistic SL SEE we have 
numerically evaluated the quantity given by Eq. (|1^) for the non-relativistic and relativistic 
regime. This approach will also provide a measure of the accuracy of the non-relativistic 
SL SEE analytic formula given by Eq. ([TSD . 

In Eig. 5 we have plotted the sudden limit SEFs for the most important astrophysical 
nuclear reactions of the CNO bi-cycle with respect to the center-of-mass energy E in the 
region of the solar Gamow peaks [Egp) with and without relativistic corrections (middle 
and lowermost curves, respectively), obtained by integrating numerically the penetration 
factors appearing in Eq. (|l^). In the same plot we have also included the non-relativistic 
SEE (uppermost curves) obtained through the analytic formula given by Eq. (|1^). We 
observe that the analytic formula overestimates the (CNO) SEFs by roughly 2% at the 
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Gamow peak energies, with the discrepancy becoming larger at lower energies. In nu- 
clear astrophysics experiments, as far as astrophysical factors are concerned, energies lower 
than the Gamow peaks are not needed for most practical purposes, unless of course there 
are suspicions about the existence of resonances. If resonances can be ruled out, or the 
experimental energies are higher than the Gamow peaks (as expected to happen in future 
experiments) then the analytic (SL) SEF formula of Eq. (|TB|) is a very reliable lower bound 
of our non-relativistic constraint. 

In the same figure, regarding the lowermost and middle curves of each reaction, although 
it is obvious that the SL relativistic curves are almost indistinguishable from the non- 
relativistic ones , their divergence is increased when moderately heavy atoms are considered 
at very small energies. Admittedly, such small effects can be neglected at the currently 
attained experimental energies (or at the even lower energies expected to be attained in the 
future). Consequently, in any case, the non-relativistic SL SEF is lower than the relativistic 
one and therefore can safely play the role of the lower bound of our constraint. 

The upper (adiabatic) limit (AL) is expected to be more reliable as it has been calculated 
under well established models. In fact the total (neutral atom) TF energy can be given by: 

oo 

E'j.^ {Z) = Z''!^ (cn + K In Z) (37) 

n=0 



where the coefficients c^, kn have been tabulated in Ref. |23|. Following the treatment used 



in Ref. |]^ the relativistic AL screening energy for proton induced reactions will be: 

oo oo 

= (Zi + 1)'/=^ [cn + fcn In (Zi + 1)] (Zi + 1)'"/=^ - ZT Y [cn + K In (Zi)] ^i"/' 

n=0 n=0 

(38 



The corresponding relativistic AL SEF can be found using Eq. ([Tq) and the above 
screening energy. In Fig. 6 we have plotted the adiabatic limit SEFs for the most important 
astrophysical nuclear reactions of the CNO bi-cycle with respect to the center-of-mass 
energy E in the region of the solar Gamow peaks [Egp) with and without relativistic 
corrections (neutral targets). Upper and lower curves always indicate relativistic and non- 
relativistic regimes, respectively. As expected the divergence between the relativistic and 
the non-relativistic curves is an increasing function of the atomic number and a decreasing 
function of the relative collision energy. The screening energy given by Eq. (^) and the 
SEFs originating from it (through Eq. (|TBp) are the highest possible values and constitute 
the upper bounds of the constraints derived in this paper 



IV. CONCLUSIONS 

The astrophysical factors S (E) for nuclear reactions ignited during advanced stages of 
stellar evolution (solar center, upper main sequence, classical novae etc.) have been ob- 
tained by performing measurements in the laboratory well above the Gamow peak, and 
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extrapolating to lower energies. Any attempt to increase the accuracy of S {E) by low- 
ering the experimental energies introduces a significant error as multi-electron screening 
effects enhance the cross sections at astrophysical energies, thus yielding an overestimated 
zero-energy astrophysical factor S (0) . The tool for correcting the respective data is the 
SEF which, in the present work, is derived both analytically and numerically for the most 
important nuclear reactions of the CNO bi-cycle by means of Thomas-Fermi model. By 
means of that model lower (SL) and upper {AL) limits are derived for the possible val- 
ues of the screening energies and the respective SEFs imposing a very reliable constraint 
on the corrections which are necessary before fitting the bare-nucleus astrophysical factor 
polynomial to the experimental data. 

Moreover, we have also investigated the effects of the degree of ionization (g) of the 
reacting atoms deriving analytically the functional dependence of SEF on it. As expected, 
the SEF is a decreasing function of q , which is again well constrained by lower and upper 
limits. Other effects studied here, are exchange, relativistic and thermal ones. The actual 
contribution of exchange effects is analysed while relativistic and thermal effects are shown 
to be negligible for most practical purposes. Finally the importance of shell effects, an 
aspect the TF model cannot take into account, is underlined by showing that the lower 
shell are mostly responsible for the screening effect while higher ones play a minor role even 
in multi-electron atomic systems. 
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FIGURE CAPTIONS 

Figure 1. The initial slope S (q) = (p (0) of the TF function (f) (x) plotted with respect 
to the degree of ionization q = 1 — N/Z of the TF positive ion. 

Figure 2. The ionic radius xq (q) of the TF positive ion plotted with respect to the 
degree of ionization g = 1 — N/ Z . 

Figure 3. The SL and AL screening energy Uq^p'^^^ for the most important astrophys- 
ical nuclear reactions of the CNO bi-cycle plotted with respect to the degree of ionization 
q = 1 — N/Z of the TF positive ion. The upper (lower) solid curve represents the en- 
hancement of the C^^ {p, 7) N^^ reaction as calculated by the TF sudden (adiabatic) limit. 
Likewise, the dashed curves stand for the N^^ (p, 7) O^^ reaction, while the dotted ones for 
the O^^ {,P,l) F^'^ reaction. Upper and lower curves always indicate sudden and adiabatic 
limits respectively. 

Figure 4. 
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The SEFs for the most important astrophysical nuclear reactions of the CNO bi-cycle 
A^^^ {p, 7) O^^ reaction with respect to the center-of-mass energy E in the region of the solar 
Gamow peaks {Egp = 27.2 keV) for various degrees of ionization q = 1 — N/Z. The lower 
(upper) solid curve represents the enhancement of the reaction as calculated by the TF 
sudden (adiabatic) limit for neutral atoms (Eqs. ([I8| ) and ^0]). Likewise, the dashed curves 
stand for a degree of ionization q = 0.3, while the dotted ones for q = 0.8. Lower and upper 
curves always indicate sudden and adiabatic limits, respectively. 

Figure 5. 

The sudden-limit SEFs for the most important astrophysical nuclear reactions of the 
CNO bi-cycle with respect to the center-of-mass energy E in the region of the solar Gamow 
peaks (Egp) with and without relativistic corrections (neutral targets). The solid curves 
represent the enhancement of the C^^ {p, 7) N^'^ reaction {Egp = 24.5 keV) with the lower- 
most and middle curves standing for the non-relativistic and relativistic SEFs, respectively, 
obtained by integrating numerically the penetration factors appearing in Eq. ([I5|) . The 
uppermost solid curve represents the non-relativistic SEE obtained through the analytic for- 
mula given by Eq. (|1^). Likewise, the dashed curves stand for the N^'^ (p, 7) O^^ reaction 
{Egp = 27.2 keV), while the dotted ones for the O^^ (p, 7) F^^ reaction {Egp = 29.8 keV) , 
with their lowermost, middle and uppermost curves having the same meaning as in the 
C^^ {p, 7) A^^^ reaction. 

Figure 6. 

The adiabatic limit SEFs for the most important astrophysical nuclear reactions of 
the CNO bi-cycle with respect to the center-of-mass energy E in the region of the so- 
lar Gamow peaks {Egp) with and without relativistic corrections (neutral targets). The 
upper (lower) solid curve represents the relativistic (non-relativistic) enhancement of the 
C^^ {p, 7) A^^"^ reaction as calculated by the TF model. Likewise, the dashed curves stand 
for the A^^^ {p, 7) O^^ reaction, while the dotted ones for the O^^ {p, 7) F^"^ reaction. Upper 
and lower curves always indicate relativistic and non-relativistic regimes, respectively. 
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Degree of Ionization (q=1-N/Z) 




Degree of Ionization: q=1-N/Z 
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